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ABSTRACT 

Angular power spectrum of the cosmic microwave background (CMB) tem- 
perature anisotropies is one of the most important on characteristics of the Uni- 
verse such as its geometry and total density. Using flat sky approximation and 
Fourier analysis, we estimate the angular power spectrum from an ensemble of 
least foreground-contaminated square patches from WMAP W and V frequency 
band map. This method circumvents the issue of foreground cleaning and that 
of breaking orthogonality in spherical harmonic analysis due to masking out the 
bright Galactic plane region, thereby rendering a direct measurement of the an- 
gular power spectrum. We test and confirm Gaussian statistical characteristic of 
the selected patches, from which the first and second acoustic peak of the power 
spectrum are reproduced, and the third peak is clearly visible albeit with some 
noise residual at the tail. 



Subject headings: cosmology: cosmic microwave background — cosmology: ob- 
servations — methods: data analysis 
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1. Introduction 

The angular power spectrum of the cosmic microwave background (CMB) temperature 
anisotropics contains a wealth of information about the properties of our Universe. The 
physics behind the shape of the power spectrum at different angular scales is well understood 
(e.g. see Hu et al.) and it therefore allows us to distinguish different cosmological models. 
The power spectrum possesses specific features, known as acoustic peaks, characterizing 
compression and rarefaction of the photon-baryon fluid around the decoupling epoch. NASA 
Wilkinson Microwave Anisotropy Probe (WMAP ) (Bennett et al. 2003a; Spergel et al. 2007; 
Hinshaw et al. 2009; Jarosik et al. 2011) has produced results that has ushered in the era of 
"Precision Cosmology", including the angular power spectrum (Hinshaw et al. 2003, 2007; 
Norta et al. 2009; Larson et al. 2011), from which cosmological parameters are estimated to 
a high precision (Spergel et al. 2003, 2007; Komatsu et al. 2009, 2011). 

In order to retrieve the angular power spectrum, however, one has to separate the fore- 
ground contamination from our own Galaxy and extra-galactic point sources in the observed 
data (Bennett et al. 2003b; Hinshaw et al. 2007; Gold et al. 2009, 2011). The standard treat- 
ment of eliminating Galactic diffuse foreground is via multifrequency cleaning, via minimum 
variance optimization for extracting the angular power spectrum, or through known fore- 
ground templates. Another issue arising from the foreground contamination is the strong 
emission of the Galactic plane, for which various masks are adopted by WMAP science 
team. Masking procedure and incomplete sky coverage thus breaks the orthogonality in the 
spherical harmonic analysis, which requires additional attention in obtaining the spherical 
harmonic coefficients (Hivon et al. 2002; Oh, Spergel & Hinshaw 1999; Ansari & Magneville 
2010). 

Apart from WMAP science team, only a couple of papers are devoted to extraction 
of the CMB power spectrum from raw data (Saha, Jain & Souradeep 2005; Saha, Prunet, 
Jain & Souradeep 2008; Samal et al. 2010; Basak & Delabrouille 2011). They all adopt the 
methodology of internal linear combination and implement quadratic minimization, which 
not only minimizes the foreground contamination, but also subtracts the power that is related 
with the chance correlation between the CMB and foregrounds (Chiang, Naselsky & Coles 
2009). 

In this paper we present a simple method of direct measurement of the CMB angular 
power spectrum from WMAP raw data, the frequency band maps. By "direct" we mean 
circumventing any foreground subtraction techniques and avoiding the issue of incomplete 
sky coverage. We also use WMAP frequency band maps, which are made possible for 
power spectrum extraction after Chiang & Chen (2011) estimate the corresponding window 
functions. 
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This paper is arranged as follows. In Section 2 we review the flat sky approximation and 
then we discuss the issue of foregrounds, instrument noise and window function in Section 
3. We test the Gaussianity of the patches taken from WMAP data in Section 4. We then 
employ the method in Section 5, and the Discussion is in Section 6. 



2. Flat sky approximation 

Standard treatment for whole-sky CMB spectral analysis is via writing the temperature 
anisotropics as a sum of spherical harmonics Y em : T(9,<p) = E«I!m fl ^m(^)' where 
9, (p are the polar and azimuthal angle, a^ m = \a,e m \ exp(i(/>^ m ) is the spherical harmonic 
coefficient and <fii m is the phase. The strict definition of an isotropic GRF requires the 
real and imaginary part of the ae m mutually independent and both Gaussian, but a more 
convenient definition is that the phases are uniformly random on the interval [0, 2ir]. The 
power spectrum can be estimated Q = (2£ + l)^ 1 Ylm \ a £m\ 2 - 

To estimate the power spectrum from small square patches, however, one can use Fast 
Fourier Transform (FFT): 

T(r) = J2 a k ex P 
k 

where r = (9, tp) and k = (kg, k v ) if the patches are chosen on the equator with sides aligned 
with the spherical coordinates. The power spectrum from the patch is = (|a^,| 2 ), where 
the angle brackets denote average for integer k over all \ayj\ 2 for k — 1/2 < \k\ < k + 1/2. 
The scaling relation between Fourier wavenumber k and multipole number £ is t = 2nk/L, 
where L is the patch size. The angular power spectrum Cf at multipole number £ is scaled 
from Cfc at Fourier wavenumber k via 

Ci = 2irk/L = L 2 Ck- (2) 

The scaling relation can be easily understood as follows if the signal is white noise. For 
spherical harmonic analysis C e = ATra^ ky /N sky whereas for FFT on a patch taken from the sky 
°patch = Cfc-Wpatch, where TVgky is the total pixel number of the sphere and A^ patc h pixel number 
of the patch. Since white noise is homogeneous, o"p atch = of ky , then Ci = 4:iiCkN patc h/ 7V B ky = 
CfcL 2 . Note that the scaling relation can be applied with minimum error for patches centred 
at 9 = 7r/2 if one uses non-equal area pixelization scheme. 

According to the scaling relation, the largest scale (smallest £) at which one can obtain 



2m(r ■ k) 
N 



(1) 
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Fig. 1. — Test of the scaling relation Eq.(2) in flat sky approximation. We simulate a full sky 
CMB map with WMAP best-fit ACDM model and take 50 patches of 24° x 24° square with 
pixel size 3 arcmin. The power spectra of the patches are scaled with a factor (27r/15) 2 and 
the sampling interval of the multipole numbers A£ = 15. The dots are from 100 patches and 
the mean is denoted by big blue dots. One can see the mean power spectrum after scaling 
fits nicely with the input one and the issue of discontinuous boundary condition is negligible. 
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the power is £ m i n = 2-k/L 1 , then the multiple numbers are sampled with the interval A£ = 
2ti/L down to the smallest scale decided by the size of the pixel p : £ max = ir/p. So the 
disadvantage of estimation of power spectrum from patches is that sampling interval A£ is 
much larger than 1, which can be viewed as intrinsic binning. 

In Figure 1 we test the scaling relation of Eq.(2). We simulate a full sky CMB map with 
WMAP best-fit ACDM model and take 50 patches, each 24° x 24° with pixel size 3 arcmin. 
One can see that the mean power spectrum from the 50 patches fits nicely with the input 
power spectrum and the error from the discontinuous boundary condition usually present in 
data analysis of square patches is negligible. 



The signal T v in the sky at frequency v is a combination of the CMB signal T c and dif- 
fuse foregrounds (synchrotron, free-free and dust emission) plus extragalactic point sources, 
altogether denoted as total foreground F u . They are measured with an antenna beam B u : 



where <g> denotes convolution and N v is the instrument noise. In order to reach the CMB 
power spectrum, we discuss below the 3 parts in Eq.(3): foreground contamination, noise 
and the window function. 



NASA Cosmic Background Explorer has measured with 10° FWHM the CMB tempera- 
ture fluctuation at a level 10~ 5 (Smoot et al. 1992). From Eq.(3) the variance of the measured 
T v includes foreground component: a 2 u = a 2 c + a% v + Cov[T c sm , T™ 1 ] + Cov[T c sm , N v ] + 
Cov[T® m , N v ], where a 2 , a\ v and a\ are the variance of the beam-convolved CMB, beam- 
convolved foreground and noise at frequency u, respectively, and the last 3 terms denote 
their covariances. For an ensemble of small patches, the average 



3. Directly Retrieving the CMB power spectrum 



T u = (T c + F u ) ®B U + N t 



(3) 



3.1. Foreground contamination and dispersion threshold 



(<$) = + (O + <4„> > (<>>) + (O- 



(4) 



1 Usually one associates multipole number t to a characteristic angular scale w on the sphere via £ = ir/vv 
because a characteristic angular scale (e.g. a scale between a cold and a hot spot) is half of one full wavelength, 
i.e. w = L/2. 




Fig. 2. — The histogram of CMB fluctuation from 1000 24° x 24° patches, which are taken 
from simulated full-sky maps with beam FWHM 19 arcmin. The mean is at 88.34/iK. 



-7- 



The CMB fluctuations (and noise) always persist in each patch, but those from the foreground 
do not. Thus we can choose patches with lower variances, as they contain less foreground 
contamination, thereby providing better estimation of the CMB power spectrum. One can 
therefore use dispersion threshold a th for controlling foreground contamination level among 
patches. 

Power spectrum is a spread of the variance into different scales, so from Eq.(4) (C£) = 
((Cfe) + + (f^), where C£ is the power spectrum of the CMB, is the window 

function, C£, and n u k are the power spectrum of the band signal, total foreground and 
instrument noise of frequency v at wavenumber k = £L/2n, respectively. For WMAP V and 
W band where the CMB dominates over the foreground outside the Galactic plane, patches 
with variance below a threshold al < a^ h shall give 

{CZ)^(Ct)WZ + {nl). (5) 

We simulate WMAP V band (i.e. with beam FWHM 21 arcmin) full-sky CMB maps 
and take in total 1000 patches of 24° x 24° square and plot in Fig. 2 the histogram of the 
dispersion a. The mean lies at 88.34/iK, which provides an indication of our choice of the 
dispersion threshold. 

There is concern that by choosing low-variance patches we are choosing CMB quiet 
area, which might result in lower power spectrum. In Fig. 3 we demonstrate that unless a 
significantly low threshold is chosen (which would result in few patches), using low variance 
as a criterion still provides a fair sample for power spectrum estimation. 



3.2. Cross-power spectrum to eliminate noise 

To eliminate the noise after choosing patches with low variance, we can employ cross- 
power spectrum on the same patch of the sky at different frequency bands. Cross-power 
spectrum (XPS) is a quadratic estimator between two maps (or patches) a and b, whose 
Fourier modes are and by,: 

< = l(( a % b k + b k a k))' (6) 

where * denotes complex conjugate and the angle brackets have the same notation as in 
Eq.(2). The advantage of XPS as an unbiased quadratic estimator for power spectrum 
estimation lies in the fact that XPS returns with its usual power spectrum (\a^\ 2 ) if a and b 
are of the same signal. If a and b are uncorrected then XPS reduces the signal by (Chiang 
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& Chen 2011) 

K _!_, (7) 

\/ A k B k \ / 27rk 

where A k and B k is the power spectrum of signal a and b respectively. The decreasing of 
the uncorrected signal is inversely proportional to the square root of the number of random 
walk, and it can be further decreased by l/y/LN with binning L = Al multipole numbers 
and averaging from N sets of XPS. Therefore XPS is useful in reducing uncorrelated signals 
while preserving the correlated one, which is employed by WMAP to extract CMB spectrum 
by crossing the foreground-cleaned maps from Differencing Assemblies (DA) (Hinshaw et al. 
2003, 2007; Norta et al. 2009; Larson et al. 2011). 

For patches on V and W band map with low variances, hence satisfied Eq.(5), we can 
write aY = a c ,bV + n\ and af = a\K + n?, where a c i is the Fourier mode of CMB, b v , 
and bj^ are that of V and W band beam, and and nj^ that of V and W band noise, 
respectively. In XPS the correlated signal (l a jL| 2 &)L&ju) * s wn& t we look for whereas those 
uncorrelated terms between CMB and noises X k cw , X k cv and between noises X k vw shall be 
decreased according to Eq.(7). 



3.3. Window Functions of the Frequency Band Maps 

The window functions of the WMAP DA maps are directly measured from Jupiter 
(Page et al. 2003; Hill et al. 2009) and are available at the official website 2 . The fre- 
quency band maps, however, are combined from the DA maps, so the corresponding window 
functions do not exist. Note that the window functions of the DA maps even at the same 
frequency band have different profiles, particularly for the W frequency band. It is then 
demonstrated in Chiang & Chen (2011) that the window functions of the frequency band 
maps can be estimated from bright point sources and it is shown that the window function 
of the W band map takes the form of that of Wl DA, whereas that of V band takes that of 
VI or V2 DA. 



4. Gaussianity of the patches 



The simplest inflation theory predicts the CMB anisotropies, amplified from quantum 
fluctuations, constitute a Gaussian random field (GRF) (Bardeen et al. 1986; Bond & Ef- 



2 http://www. lambda. gsfc.nasa.gov/product/map/current/ 



-9 - 



stathiou 1987). If the CMB is indeed statistically isotropic Gaussian, the angular power 
spectrum furnishes a complete statistical description. One of the properties of GRF is its 
phases from harmonic analysis are uniformly random in [0, 2ir]. Based on the random phase 
hypothesis we can test Gaussianity of the selected patches by employing the Shannon en- 
tropy of Fourier phases S = — ^2pi lnpi5<f), where PiS(j) is the distribution probability at the 
ith interval in [0, 2tt] and J2piS(p = 1. It can be used to test for uniformity: p = p{4>k) and 
independence (non-association): p = p(D) where D(Ak) = (fik+Ak — <t>k (Chiang & Coles 
2000; Coles & Chiang 2000). If the phases are uniformly random, then p^ = (27r) _1 and the 
entropy reaches a maximum S max = \n2n. 

5. Angular power spectrum of the CMB from the WMAP frequency band 

maps 

In this section we apply our method on WMAP frequency V and W band maps to extract 
the CMB angular power spectrum via Fourier analysis on 24° x 24° patches. Although one 
can take patches with a smaller size rendering more patches from the whole sky, it would 
increase the noise power spectrum level, and consequently, the XPS residual shown in Eq.(5). 

We first take WMAP V band and choose 47 patches with a < 98/iK (after deleting 
bright point sources exceeding 5a of the patch). Note that in Fig. 2 the mean of the 1000 
patches is a = 88.34/iK, but those taken in real maps with pixel noise have higher values. 
Before extracting the power spectrum, we test the Gaussianity of the 47 patches by employing 
Shannon entropy of Fourier phases for their association. In Fig. 4 we show the normalized 
histogram of Shannon entropy for association of phases for AA; = 1 on the left and 2 on the 
right. We also plot the histogram from 1000 24° x 24° Gaussian patches. 

We then calculate the XPS from the same patches of the V and W band to eliminate 
the noise. They are then deconvolved by (\/WyWw)~ 1 , where Wy and Wy? are the window 
function of the V and W frequency band map, respectively. We show in Fig. 5 the retrieved 
angular power spectrum. One can see our simple method yields the CMB power spectrum 
with clear 1st and 2nd Doppler peak, which match the result from WMAP science team. 
The 3rd peak is also visible albeit with higher amplitude at the tail than theirs, which is due 
to the residual from XPS. 
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Fig. 3. — Demonstration of CMB quiet areas from low variance. From a simulated V band 
CMB map (high I are smoothed by beam FWHM 21 arcmin), we set different threshold cr th 
to see if choosing CMB quiet areas affects the power spectrum estimation. On the left panel 
we choose threshold a th = 82/iK with only 3 patches in the plot. The mean power spectrum 
(big blue dot) is indeed lower than the input one (solid line), particularly for low £. In the 
middle panel we plot the 9 patches that have a lower than 84/iK. On the right panel with 
cx th = 88/iK (the mean from Fig. 2), one can see the mean power spectrum from 34 patches 
fits well with the input one. 
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Shannon Entropy of phase association for Ak=1 
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Fig. 4. — Normalized histogram of the Shannon entropy for Fourier phase association from 
the 47 selected patches. The phases of the Fourier modes equivalent to multipole number 
£ < 1050 are taken for calculation. The Shannon entropy S of phase association between 
Ak = 1 is shown on the left panel and Ak = 2 on the right, the normalized histogram for 
the 47 patches is shown in solid curves. For comparison we plot in dash curve the histogram 
from 1000 24° x 24° patches taken from full-sky Gaussian maps. 
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Fig. 5. — Direct measurement of the CMB angular power spectrum. From WMAP V band 
map we choose patches with a < 98/iK (after eliminating bright point sources), and we take 
the cross-power spectra of patches between WMAP V and W band. After deconvolution of 
the window functions, the power spectra of the 47 patches are shown in black dot and the 
mean power spectrum is in big blue dot. For comparison we plot in big orange dot the power 
spectrum binned (A£ = 15) from that by the WMAP science team. The best-fit ACDM 
model is in solid line. 



6. Discussion 



Since WMAP data release, full-sky analysis has become the standard way for power 
spectrum estimation as it has offered for the first time detailed measurement at super-horizon 
scales. The method we present in this paper utilizes small patches, hence has intrinsic 
limitation on the largest scale we can measure. It nevertheless adopts a totally different 
methodology and provides a more intuitive way to obtain the power spectrum on all but the 
very large scales. This method can be readily applied on the upcoming Planck data. 

We acknowledge the use of Healpix 3 package (Gorski et al. 2005) and the use of Glesp 
4 package. The author would like to thank Peter Coles, Andrew Jaffe, Dipak Munshi and 
Keiichi Umetsu for useful discussions and suggestions. 
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